Exact and quasi-resonances in discrete water-wave turbulence 
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The structure of discrete resonances in water-wave turbulence is studied. It is shown that the 
number of exact 4-wave resonances is huge (hundreds million) even in comparatively small spectral 
domain when both scale and angle energy transport is taken into account. It is also shown that angle 
transport can contribute inexplicitly to scale transport. Restrictions for quasi-resonances to start 
are written out. The general approach can be applied directly to mesoscopic systems encountered 
in condensed matter (quantum dots), medical physics, etc. 
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1. Introduction. Statistical wave turbulence the- 
ory deals with the fields of dispersively interacting waves. 
Examples of these wave systems can be found in oceans, 
atmospheres, plasma, etc. Interactions between waves 
are similar to interactions between particles and can be 
described by kinetic equations (KEs) analogous to KE 
known in quantum mechanics since 1930th. Wave KE is 
in fact one limiting case of the quantum Bose-Einstein 
equation while the Boltzman kinetic equation is its other 
limit. First wave KE for surface gravity waves has been 
presented in [l[ while general method of construction KEs 
for many other types of waves can be found in Q ■ One 
of the major achievements of the statistical approach is 
establishing of the fact that power energy spectra (simi- 
lar to famous Kolmogorov 5/3 law) are exact stationary 
solutions of kinetic equations 0] . The limitations of sta- 
tistical turbulence theory are due to the fact that it does 
not describe spatial unevenness of turbulence, i.e. orga- 
nized structures extending over many scales, like boulders 
in a waterfall, remain unexplained. Appearance of these 
structures is attributed to mesoscopic regimes which are 
at the frontier between classical (single waves/particles) 
and statistical (infinite number of waves/particles) de- 
scription of physical systems. Mesoscopic (effectively 
zero-dimensioned) systems is very popular topic in var- 
ious areas of modern physics - from wave turbulence to 
condensed matter (quantum dots, to medical physics. 
For instance, in [f| dynamics of blood flow in humans is 
studied, cardiovascular system is described by a few cou- 
pled oscillators and synchronization conditions are inves- 
tigated. Synchronization or resonance conditions have 
the same general form for wave and quantum systems 
(see, for instance, 0] for 4-photon processes); and have 
to be studied in integers. Just for simplicity of presen- 
tation we prefer to stay with wave terminology while all 
examples in this Letter are taken from wave turbulent 
systems. 

Resonance conditions have the form 

0i±0 2 ±...±0 s = 0, £i±fc 2 ±....±fe s = o (1) 

where 0j = 0(fcj), s < 00, with fc and 0j = 0(fcj) 
being correspondingly wave vector and dispersion func- 
tion. Specific features of these systems described by 



Fourier harmonics with integer mode numbers were first 
presented in [?J (we call them further discrete wave sys- 
tems, DWS). A counter part to kinetic equation in DWS 
is a set of few independent dynamical systems of ODEs 
on the amplitudes of interacting waves. Mathematical 
theory of DWS was developed in Q with general under- 
standing that discrete effects are only important in some 
bounded part of spectral space, |fc| < ko, with some 
small finite fco , while the case |fc| >> ko is covered 
by kinetic equations and power-law energy spectra. This 
general opinion was broken recently as result of numerical 
simulations with Euler equations for capillary waves Q 
and for surface gravity waves [Io| where discrete clusters 
of waves were observed simultaneously with statistical 
regime. Moreover, experimental results ll| show that 
discrete effects are major and statistical wave turbulence 
predictions are never achieved: with increasing wave in- 
tensity the nonlinearity becomes strong before the system 
loses sensitivity to the fc-space discreteness. 

2. Discrete wave systems. In [13] a model is pre- 
sented which explains the appearance of discrete wave 
clusters in the large spectral domains with |fc| >> fco. 
It shows that energy power spectra are valid not in all 
spectral domain with big fc but have " holes" all over the 
spectrum, in some integer fc, which describe discrete dy- 
namics of mesoscopic regimes. This understanding put 
ahead a novel computational problem - computing inte- 
ger solutions of |T|) in the large spectral domains. Indeed, 
((T|) turns into 

\fk~i + \fk~2 = \fk~3 + \fk~i , k 1 +k 2 = k 3 + k 4 , (2) 

for 4-wave interactions of 2D-gravity water waves, where 
fc, = (nii,ni), Vi = 1,2,3,4 and ki = \kj\ = \fm% + ri\. 
This means that in a domain, say |m|, \n\ < D ~ 1000 , 
direct approach leads to necessity to perform extensive 
(computational complexity D 8 ) computations with inte- 
gers of the order of 10 12 . The full search for multivariate 
problems in integers consumes exponentially more time 
with each variable and size of the domain to be explored. 
The use of special form of resonance conditions allowed 
us to develop fast generic algorithms [l3| for solving sys- 
tems of the form ([I]) in integers. The main idea of the 
algorithms for irrational dispersion function is that if vec- 
tors ki construct an integer solution of fl}, then at least 
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for some io,jo the ratios 0i o /0j o have to be rational num- 
bers. Some other number-theoretical considerations were 
used in case of rational dispersion functions. In particu- 
lar, all integer solutions of ^ in domain \m\, \n\ < 1000 
were found in a few minutes at a Pentium-3 (cf.: direct 
search for them in smaller domain \m\, \n\ < 128 took 3 
days with Pentium-4 

Using our programs we studied how the structure of 
discrete resonances depends on the form of dispersion 
function and on the chosen s. The interesting fact is 
that characteristic structure is the same for different dis- 
persion functions (examples with = m/n(n + 1), = 
1 / y/m 2 + n 2 , = m/(n 2 + m 2 + 1 ) and others were stud- 
ied) if s = 3. Most waves do not take part in resonances 
and interacting waves form small independent clusters, 
with no energy flow between clusters due to exact reso- 
nances. Our conclusion about wave systems with 3-wave 
interactions is therefore that exact resonances are rare 
and quasi-resonances, i.e. those satisfying 

0i ± 02 ± ... ± 0s = n > 0, k 1 ±k 2 ± .... ± k s = (3) 

with s = 3 can be of importance for some applica- 
tions. The situation changes substantially in the case 
s > 4 which is illustrated below with surface gravity 
waves taken as our main example. 

3. Exact resonances. The major difference be- 
tween 3- and 4-wave interactions can be briefly formu- 
lated as follows. Any 3-wave resonance generates new 
wave lengths and therefore takes part in energy trans- 
fer over scales. In a system with 4-wave resonances 
two mechanisms of energy transport are possible: 1) 
over scales, if at least one new wave length is gener- 
ated, and 2) over angles, if no new wave lengths are 
generated. Examples of these two types of solutions 
for (0) are (-80, -76) (980, 931) =*> (180, 171)(720, 684) 
and (-l,4)(2,-5) =*> (-4, 1)(5, -2), we call these so- 
lutions further scale- and angle-resonances correspond- 
ingly. Two mechanisms of energy transport provide sub- 
stantially richer structure of resonances and the num- 
ber of exact resonances grows enormously when com- 
pared to 3-wave resonance system. Say, for dispersion 
= l/\Zm 2 + n 2 and 3-wave interactions there exist 
only 28156 exact resonances with \m\, \n\ < 1000 while 
for dispersion = (m 2 + n 2 ) 1 / 4 and 4-wave interac- 
tions in the same domain the overall number of exact 
resonances is about 600 million. However, major part 
of these resonances are angle-resonances. In the do- 
main |m|,|n| < 1000 we have found only 3945 scale- 
resonances, i.e. transport over the scales is similar to 
the case of 3-wave interactions. Some isolated quar- 
tets do exist, also among angle-resonances, for instance 
(— 1, — 1)(1, 1) (—1, 1)(1, —1) but they are rather rare. 
Very important fact is that angle and scale energy trans- 
port are not independent in following sense. One wave, 
say (64,0), takes part in 2 scale- resonances one of them 
being (64,0) (135, 180) =S> (80, 60)(119, 120) and the 



wave (119,120) takes part then in 12 angle-resonances, 
and further on (see Figs Jll2|) . 

It is important to understand that for s > 4 1) no 
new type of resonances appear, and 2) existence of angle- 
resonances depends on the placing of signs in the first 
equation of (TTJ): for instance, if it has the form 0i = 
02 + ... + 0s, no angle-resonances are possible. 




Figure 1: Wave (64,0) takes part in 2 scale-resonances. The 
upper number in the circle is m and the lower - n, and red 
thick lines drawn between vectors on the same side of the 
Eqs.©. 



12 symmetric solutions 




Figure 2: Wave (119,120) takes part in 12 angle-resonances. 

The natural question now is whether quasi-resonances 
are in fact of importance in a wave system possessing 
such enormous number of exact resonances. 

4. Quasi-resonances. From now on we are inter- 
ested in discrete quasi-resonances which arc integer so- 
lutions of (J3j> with some non-zero resonance width fi. 
Notice that co : Z x Z — > Q(<?) , where q is an al- 
gebraic number of degree 4 and the field Q(q) de- 
notes corresponding algebraic expansion of Q . To es- 
timate a linear combination of different 0^ over Q we 
use generalization of the Thuc-Siegcl-Roth theorem [l5j |: 
If the algebraic numbers ai, a.2, a s are linearly in- 
dependent with 1 over Q , then for any e > we 
have \piOt\ + P2«2 + •■• + Ps<^ s — p\ > cp~ s ~ e for all 
p,Pi,P2, ...,Ps G ^ with p = maxi|p 2 ;|. The non-zero 
constant c has to be constructed for every specific set 
of algebraic numbers separately. For = (to 2 +7?, 2 ) 1 / 4 
and four-term combination this statement means in par- 
ticular that 1 0i ± 02 ± 03 ± 0a\ > 1 if at least one of 
0i is not a rational number. As it was pointed out in 
[131 ] all integer solutions of ((2]) have one of two forms: I) 
0i = W f \ V* = 1, 2, 3, 4 or II) 0i = 3 = 7K71 74 , 2 = 
04 = 72<72^ 4 ' 9i 92- Here ^i,q, qi,Q2 are some integers 
and g, (71,(72 have form p^p^ 2 ■■■Pn n > with all different 
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primes pi,...p n while the powers e\, ...e„ € NU {0} arc 
integers all smaller than 4. It follows that for all wave 
vectors but those of the form I) with q = 1, there exists 
a global low boundary for resonance width, fi > 1, 
necessary to start quasi-resonances. In the spectral do- 
main |m|,|n| < 1000 only 136 scale-resonances do not 
have global low boundary. But even for them local low 
boundary exists - defined by the spectral domain T = 
{(m,n) : < |m|,|n| < D < oo} chosen for numerical 
simulations. Indeed, let us define Qjj = min p where 
n p = \ (k p 1 )±0(k p 2 )±...±0(k p 4 )\, k p = (roj.nf) e T, for 

Vj = 1,2,3,4, 0(k p )±0(k p )± ...±0(k p ) ^0 Vp, and 
index p runs over all wave vectors in T , i.e. p < AD 2 . 
So defined fl p obviously is a non-zero number as a min- 
imum of finite number of non-zero numbers and Qd 
is minimal resonance width which allows discrete quasi- 
resonances to start, for chosen D. 

Step of numerical schema < S <C 1 is another 
parameter important for understanding quasi-resonant 
regimes, and inter-relation between S, fljj and f2 
describes them all. For instance, if fi^ > f2, any chosen 
5q > will allow some number of quasi-resonances, 
sa y N (S >n D ) , and for any 6 > S => N^>n n ) ^ 
N(s >n D )- On the other side, if 5 is decreasing to 0_d 
from above, 5 — > iljj+O, the number of quasi-resonances 



reaches some constant level N„ 
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N (s=n D ) 



If S is increasing to from below, (5^0^ — 0, the 
number of quasi-resonances is N m i n = const . 

This fact has been first discovered in the numerical 
simulations [l6j |. both for capillary and surface gravity 
waves (maximal spectral domain studied was < m < 
2047, < n < 1023 and only scale-resonances were re- 
garded). In case of gravity surface waves, increasing of 
from 10~ 10 to 10~ 5 does not changes number N of quasi- 
resonances while increasing of il from 10~ 4 to 10~ 3 
yields increasing N — > lOiV. It turned out that limiting 
level N m i n is for capillary and and ~ 10 5 for surface 
waves, for the same discretization by a rectangular mesh. 
It was attributed to the following fact: quasi-resonances 
are formed in some vicinity of exact resonances which do 
exist in the case of gravity waves and are absent in the 
case of capillary waves. 

5. Topological structure of resonances. Graph- 
ical presentation of discrete 2D-wave clusters suggested 
in [8| is to regard each 2D-vector as a node of integer 
lattice in spectral space and connect those nodes which 
construct one solution (triad, quartet, etc.) We demon- 
strate it in Figl3] (upper panel) taking for simplicity 3- 
wave interactions with = 1/ \fm 2 + n 2 (ocean planetary 
waves). Obviously, geometrical structure is too nebulous 
to be useful. On the other hand, topological structure 
in Fig[3] (lower panel) is quite clear and gives us imme- 
diate information about the form dynamical equations 
covering behavior of each cluster. The number in brack- 
ets shows how many times corresponding cluster appears 



in the chosen spectral domain. All similar clusters are 
covered by similar systems of ODEs (written for simplic- 
ity for real- valued amplitudes): Ax = ctiA 2 A 3 , A 2 = 
a 2 AiA 3 , A 3 = a^A 1 A 2l in the case of a "triangle" 
group, two coupled systems of this form in the case of 
"butterfly" group and so on. A 3- wave system has been 
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Figure 3: Upper panel: Geometrical structure, D = 300. 
Lower panel: Topological structure, D = 50. The number in 
brackets shows how many times corresponding cluster appears 
in the chosen spectral domain. 

chosen here as an illustrative example for its simplicity. 
They have only one type of vertices in their graphical 
presentations: nodes, corresponding to exact resonances. 
Some graph-theoretical considerations allow to construct 
isomorphism between topological elements of the solution 
set of ([1]) and corresponding dynamical systems, for the 
case s — 3. This yields a constructive method to gen- 
erate all different dynamical systems in a given spectral 
domain. For instance, all graphs shown in Fig[3l lower 
panel, are described by 4 different dynamical systems; for 
all isomorphic graphs, corresponding dynamical system 
has the same form though coupling coefficients Oj have 
different magnitudes, of course. Some programs are writ- 
ten in MATHEMATICA which allow for a few specific 
examples to a) find all integer solutions of ([T]), b) gener- 
ate their geometrical and topological structure, c) write 
out explicitly all corresponding dynamical systems. Only 
small number of these systems are known to be solvable 
analytically (mostly those corresponding to clusters of 3 
to 5 waves only) while larger systems should be solved 
numerically, of course. Knowledge of specific form of a 
dynamical system allows in many cases to write out some 
conservation laws and thus simplify substantially further 
numerical simulations. 

In case of s-wave interactions with s > 4 construc- 
tion of a corresponding graph must be substantially re- 
fined: a graph with 3 different types of vertices should 
be constructed, corresponding to the waves taking part 
in 1) angle-resonances, 2) scale-resonances, 3) both types 
of resonances, in order to provide simultaneous isomor- 
phism of graphs and dynamical systems. This work is 
under the way. 
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Knowledge of the resonances structure might con- 
tribute to short-term forecast of wave field evolution, 
for in direct numerical simulations [3],[l3| discrete res- 
onances were observable not at the time scale 0(t/e A ) of 
kinetic theory (e is steepness of wave field) but at linear 
time scale t. 

6. Conclusions. 

1. 3-wave systems can possess only scale- resonances 
which are rare, in this case quasi-resonances might be of 
importance in energy transport; 

2. s-wavc systems with s > 4, depending on the sign 
setting in {TJ), may also posses angle-resonances which 
contribute inexplicitly into energy transport. In systems 
like ([2]) where angle-resonances are allowed, there exist 
hundreds million of exact resonances in a comparatively 
small spectral domain \m\, \n\ < 1000. 

3. For some polynomial irrational dispersion function, 
global low boundary £1 for quasi-resonances to start 
can be computed which 1) does not depend on the cho- 
sen spectral domain, and 2) is valid for the most part of 
exact resonances in a system. If dispersion is a rational 
function, only local low boundary £lrj exists (it follows 
from the fact that Q is dense everywhere in R). 

4. Any interpretation of results of numerical simula- 
tions with dispersive wave systems has to take into ac- 
count interplay between 5, Qo and fl. 

5. Specially developed graph presentation of the solu- 
tion set of fl} allows to construct isomorphism between 
independent cluster of resonantly interacting waves and 
corresponding dynamical systems. This yields construc- 
tive algorithm for generating dynamical systems symbol- 
ically, for instance using MATHEMATICA. 

6. The same approach (the use of the algorithms from 
low boundary computing, graph construction and 

generation of dynamical systems) can be used directly 
for any mesoscopic system with resonances of the form 
(fT]) or, more generally, of the form 

Pl01±p 2 02±...±P S 0s = 0, Plkl ±p 2 k 2 ± ....±p s k s = 

with integer pi , in this case global low boundary will 



depend on max^ \pi \ . 
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